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Abstract 

The influence of charged STM tip on the electron transport through quantum 
states on a surface is studied both theoretically and experimentally. The 
current and the differential conductance calculations are carried out by means 
of the Green's function technique and a tight-binding Hamiltonian. It is 
shown that sharp STM tip is extra occupied and this additional charge breaks 
the conductance symmetry for positive and negative STM voltages. The 
experiment on Ag islands with two STM tips (blunt and sharp) confirms our 
theoretical calculations. 
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1. Introduction 

The invention of scanning tunnelling microscope (STM) in the early 
1980's [H revolutionized imaging processes of various surfaces. This instru- 
ment allows us to study exact positions of atoms (in conducting surface or 
in small nanostructures at the surface) with atomic resolution. Moreover, 
the electrical properties (tunnelling spectroscopy, STS) can be investigated 
as well. 

The STM resolution depends on the tip quality: the sharper tip the better 
resolution can be achieved and only one atom at the end of the tip is the 
ideal case. The last tip atom properties are crucial for electron transport 
through the STM, 0, 0, 0, 0]. Using such a sharp tip the spectroscopy or 
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topography studies of atomic objects and quantum states on a surface are 
possible with the best resolution. In order to obtain useful information on 
density of states of investigated objects an appropriate treatment of STS 
results is required. The most popular theoretical description of the STS 
transport have originally been developed by Tersoff and Hamann 0] . Within 
the WKB approximation the tunnelling current is expressed in terms of the 
surface and tip local density of states (DOS) and the transmittance through 
the system. However, the energy convolution between both density of states 
makes it difficult to interpret STS results, e.g. [3,0]. It is expected that the 
structure of the tip density of states should be reflected in STS results and 
can significant change the position or intensity of the surface DOS peaks, 0|. 
Moreover, new maxima can appear on the differential conductance which 
are not connected with the surface DOS, e.g. [4j. It is also known that 
metal tips induce a band bending on the semiconductor surface (tip-induced 
band bending effect), e.g. [g|, 0, [Io[, which strongly influences the tunnelling 
current and leads to new peaks in STS results. 

In many STM experiments on single atoms, molecules or even flat atomic 
surfaces the STS results (differential conductance curves) obtained for pos- 
itive and negative STM voltages are asymmetrical. Also the topography 
images for both voltages are often quite different, cf. 0, EH. The main rea- 
son of such behavior is asymmetry in the density of states of investigated 
objects. Such an object can be characterized itself by asymmetrical DOS 
(e.g. highly asymmetric atomic structure in the tunnelling current between 
neighboring graphite atoms was observed in Ref. Il2j) or this asymmetry can 
appear e.g. due to the coupling object-surface, [131 ] . Moreover, non ideal tip 
geometry can also lead to different pictures for positive and negative volt- 
ages. Note, that in real STM experiments it is almost impossible to obtain 
fully symmetrical results. 

From the basic electrodynamic rules, one expects large charge concentra- 
tion on sharp or curved conducting materials, like e.g. on STM tips. Charge 
occupation of such tips is very important as concerns electrical properties 
and can lead to asymmetry in the differential conductance (versus the posi- 
tive and negative STM voltages) even for objects which are characterized by 
fully symmetrical density of states. We show in this paper that the sharper 
tip, the more charge is accumulated on it and the conductance symmetry can 
be broken in this case. (Note, that this effect is not relevant to electrostatic 
forces which appear due to differences between the Fermi levels of a sample 
and the tip or due to applied sample voltage, e.g. 10, [l3])- The goal of this 
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paper is to show that additional charge localized at the STM tip, in addition 
to other possible sources, influences the current-voltage characteristics and 
is responsible for the asymmetry effect in the STS studies. In our calcula- 
tions we consider three models of the STM tip and use the Green function 
formalism together with a tight-binding Hamiltonian to obtain the conduc- 
tance through quantum states on a surface. To illustrate experimentally our 
theoretical predictions the STS studies of thin Ag films on Si(lll)-(6x6)Au 
substrate were performed. During the measurement a sudden variation of 
the tip length took place which allowed us to compare the STS results for 
two kinds of STM tips (sharp and blunt). Note, that small changes in the tip 
structure often is observed during a series of scans, e.g. [2|. These changes 
have only a minor effect on topographic pictures but can drastically alter the 



conductance curves, cf. also [15 



The paper is organized as follows. Additional charge at the STM tip 
is analyzed in Sec. [2j In Sec. [3] the Hamiltonian for three model tips and 
general formulas for the STS current and the transmittance are shown. In 
Sec. H]the results for a single atom and short wire on a surface are discussed. 
The STM experimental results and their theoretical description are shown in 
Sec. |5j The last section, Sec. El is devoted to conclusions. 



2. Charged STM tip 

The main idea of this paper is that the STM tip is charged i.e. it possesses 
an additional charge. This charge influences the STS results and leads to 
asymmetry in the conductance or current-voltage curves. In this section we 
analyze a simple model of the STM tip and explain why an additional charge 
is localized at the last tip atom. 

First, let us analyze the position (on the energy scale) of the apex atom 
energy level, Sstm-i versus the chemical potential of the STM electrode. For 
£stm > Ustm (&STM < Ustm) the tip is not charged (is charged) and for 
£stm = Ustm is neutral. In most of our calculations we set the energy level of 
the apex atom below the Fermi energy of the STM electrode because in this 
case the STM tip is charged. Here we give the reason of our choice. Using 
the basic electrodynamic relations one can show that additional charge is 
gathered on sharp edges and curvatures. To confirm this rule for the STM 
we consider here a simple model of the STM tip i.e. one-dimensional bend 
monoatomic wire. In Fig. [T^i the schematic view of a straight wire with 
the nearest neighbor couplings between atoms, Vq, is shown. For a bend 
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Figure 1: The schematic view of (a) iV-atom straight wire coupled with the left and right 
electrodes, and (b) bend wire with nonzero hopping, V x , between x — 1 and x + 1 atoms 
in the wire. Figure (b) represents the STM tip which is fabricated by bending a straight 
one-dimensional wire (panel a). The figure depicts only a part of junction, namely the tip 
without a sample. 

wire, Fig. [lb, also the next neighbor hopping has to be considered, V x . The 
Hamiltonian for this system is given by 



N 

i=1 keL/R 

N-l 



H = ^ :,<■; r, • ^ € < rZr 



+ V k c t Cl / N + V a c t c i+i + Vxct-^x+i + h.c. , (1) 

k&L/R i=l 

where x stands for the middle atom in our wire, x = (N + l)/2, iV is odd. 
The main quantity of our interest is the occupation of this middle atom, 
n x , which represents the last tip atom. The local charge at this atom can 
be expressed by means of the retarded Green's function, G r xx i.e. n x = 
J E ^ImG r xx (s)de. The funct ion G xx can be obtained from the equation 
of motion for the retarded Green's function and in our case [for iV atom 
wire with the same electron energies, £j = Eq, and within the wide band 
approximation, T = T L = T R = 2n J2keL \^k\ 2 ^( £ ~~ £ k)] ^ becomes: 

d2 _ t/2 r>2 

QT ( {?\ — D x-\ v x JD x-2 

xxK } B N + &B N _ X - V X [2V 2 + V x (e - e )]B 2 x _ 2 1 ] 
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Figure 2: The occupation of the middle atom in the wire as a function of V x parameter 
(see Fig. [T]d) for N — 3 (solid line) and N — 5 (broken line). The other parameters are: 
r = 1, £ = 0, E F = 0, Vo = 4, T = OK (all energies in T units). 



where B x = det A° x + det and A° x matrix reads: 



( e-£ -Vo 
-V e - e -V 
-Vo £-so 



\ 



(3) 



-V e - so J v 

/ xxx 



The determinant of this matrix can be expressed by means of the Chebyshev 
polynomials of the second kind, 16|], and finally the retarded Green's func- 



tion, G r xx , can be obtained analytically. It helps us to find the occupation of 
the middle atom in the wire, as a function of V x parameter. 

In Fig. [2] we show the charge localized at the middle atom versus the 
coupling parameter V x for the wire length N = 3 (solid line) and N = 5 
(broken line). The parameter V x is responsible for the wire curvature i.e. for 
V x = the wire is straight and for nonzero V x we have a kind of bend wire 
as is depicted in Fig. [Tb. Here we consider the case e = = = which 
leads to half occupation of all wire sites for V x = (^l/r is the chemical 
potential of the L/R electrode). However, for nonzero V x the occupation of 
the middle atom increases and e.g. for V x = Vo = 4 it is about 30% larger 
than for the straight wire. This result indicates unquestionably that there 
is additional electron charge at the middle atom. This atom corresponds to 
the last tip atom and is more occupied in comparison with other atoms. 
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Here we should comment on the last result. Of course, STM tips are not 
fabricated by bending a straight monoatomic wire. However, a sharp STM 
tip, with only one atom, can be described as a kind of short bend wire e.g. 
with iV = 3 or N = 5 atoms. In this case the coupling V x is responsible 
for the tip curvature and the larger V x the sharper tip is (there is physical 
limit on V x - it cannot be larger than twice or a few values of Vq parameter). 
Thus, this approach is suitable for description of sharp tips and failed for 
rather wide ones. Note, that in the literature STM tips are often much more 



or an 



simplified and modelled by a semi-infinite chain of atoms, e.g. [171 . Il8l 
ideal fermi gas electrode (energy structure-less). 

On the other hand, one can consider the last tip atom as a kind of ad- 
ditional atom (adatom) which is put on metallic STM electrode. This atom 
cannot be treated in the same way as atoms inside the STM electrode, first 
of all because its neighborhood is different in comparison with other atoms 
in the tip. Moreover, the energy of affinity (and ionization) level at the apex 
atom depends on the distance between this atom and the STM metallic sur- 
face. This effect is known in chemisorption processes, i.e. the affinity energy 
level decreases its value versus the surface Fermi energy with the distance 
adatom-surface and is minimal for an adatom placed directly on a surface, cf. 



19j. In that case electrons from the surface can occupy the adatom affinity 
state which leads to minimalization of the total energy (cf. also the calcu- 
lated local DOS of Co atoms at Au surface which is characterized by local 
peaks under the Fermi level, 2p|). This explanation support our idea that 
the STM tip is charged. Also dynamical polarization effects (surface-adatom 
or surface-molecule) can renormalize molecular states, [2l| . 

Taking into account the above considerations we assume in our calculation 
that the energy of the last tip atom, Estm, l a Y s slightly below the STM 
chemical potential, i.e. Estm < Ustm, and thus is more occupied. One 
can obtain the same effect for Estm = Ustm but the STM tip should be 
considered as a bend wire, cf. Fig. [lb. In the last case the middle atom in 
the wire is extra occupied (its density of states possesses a local maximum 
below the Fermi energy). It allows us to treat the STM tip as a metallic 
electrode with only one atom under the condition Estm < Ustm- 

3. Theoretical STM model 

The system consists of the STM tip, an object (few-atom system, i.e. one 
atom or short wire) and a metallic surface, cf. Fig. [3] for only one atom on 
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Figure 3: Sketch of one-atom object on a surface and three STM tip models: (a) STM tip 
as a metallic electrode, (b) tip consisted of a metallic electrode with one atom (coupled 
with the STM electrode via r STM0 ); the apex atom is also coupled with the investigated 
object via to parameter, (c) the same as in (b) but for nonzero coupling between the tip 
electrode and the object via r STM1 parameter. 



the surface. The STM tip can be represented by (i) a metallic electrode i.e. 
tip 1, Fig. [3^ - electrons flow from this electrode directly to the investigated 
object, (ii) a metallic electrode with the last tip atom on it, i.e. tip 2, Fig. [3b 
- electrons flow from the electrode through this apex atom then through the 
object on the surface, and (iii) a metallic electrode with the apex atom but 
electrons can flow from the STM electrode directly to the object or through 
the last tip atom - tip 3, Fig. [3b. Tip 2 in Fig. [3] is suitable for very sharp 
STM tips and tip 1 (widely used in the literature) describes rather blunt tips. 
The third case, tip 3, can be treated as a composition of the first and the 
second tips. To obtain the current flowing from the STM electrode we use 



the Green function formalism and thigh-binding Hamiltonian, cf. [14 
which can be written as follows: 

H = H S tm + H surf + H object + V , (4) 

where H$tm describes electrons in the STM tip. For a metallic STM elec- 
trode it can be expressed in the form (ideal Fermi gas): Hstm = T]i ct cr , 

V V O J Oliu Z^tksTM KSTM ksTM k STM ' 

which corresponds to Fig. [3^ or with one tip atom, cf. Fig. [3b, c: 

H S tm = { £ k c t c k + h c t c s™ + h - c -) + £ stmCstm c stm , (5) 

fee STM 

Here ££ STM and Estm represent electron energies in the STM electrode and at 
the tip atom. c + and c are the creation and annihilation electron operators 
in an appropriate electron state and t? stands for the coupling parameter 
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between the STM electrode and the tip atom (r 5 ™ in Fig. [3] depends on 
this parameter). The surface, similarly to the STM electrode, is modelled as 
an ideal Fermi gas by the Hamiltonian, H sur t = Vr e% ct c% . The 

J — ™surf K surf fosurf sur f 

object Hamiltonian depends on the number of atoms at the surface and e.g. 
for a single atom it can be written as H bj ect = eic^ c\, but for a short wire 
it has the following form 



N N-l 



H oh ject = ^ e i C t C i + E V '< C t C i+l + h ' C - ( 6 ) 
i=l i=l 

The coupling Hamiltonian is responsible for electron tunnelling between the 
STM and the surface electrode and for H$tm, Eq. 0, and H ob j ect = eic^ci, 
(tip 2, Fig. [3b) it can be expressed as: V = J2k c t c i + to c t c STM + 

surf surf rCsurf 

h.c, where t is the hopping integral between the tip atom and the object 
atom. Note, that for iV-site wire on the surface one should change the cou- 
pling Hamiltonian and instead of c\ operator one has to sum q operators over 
all atomic sites in the wire i.e. i = 1, N. Moreover, for the system shown 
in Fig. [3b one has to add to the Hamiltonian the term which is responsible 
for the coupling between the STM electrode and the object atom. Thus the 
coupling Hamiltonian becomes: 



N 

v = E E ^/ c L/* + toc ^™ + E v * STM c L^ + hx - > ( ? ) 

i=i % . 

^surf 



where we assume that the STM electrode is connected with the first atom 
in the wire via V% TU element (it is responsible for r STM1 parameter, cf. 
Fig. [3b). The role of electron-electron interactions in the system is discussed 
in Sec. 14.21 Electron transport properties are analyzed within the framework 
of the Green's function method. The tunnelling current flowing from the 



STM electrode to the surface can be obtained from, 22 



2e 

T 



J deT{e)[f STM {e) - fsurf(e + V)] , (8) 



where f(e) is the Fermi function and the external voltage is expressed by 
means of the STM and surface chemical potentials i.e. V = ustm — Hsurf- 
The transmission function for considered here systems reads: 

T(e) = Tr{f STM G r f surf G a } , (9) 
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where TsTM/surf is the matrix composed of the coupling parameters be- 
tween the STM electrode or the surface with the object atoms. In gen- 
eral the elements of these matrixes can be written in the form (T S urf)ij — 
2vr V-t V* r 5(e — E? ), and similar for Yrtm- The next ma- 

L — JK surf l ^surf J,K sur f K surf' 

trixes in Eq. [9] are composed of the retarded and advanced Green functions 
(G a = G r *), and can be obtained from the equation of motion for the retarded 
Green's function. Note, that for our systems, Fig. [31 the transmittance can 
be obtained analytically and all matrix elements will be specified later. The 
knowledge of the current flowing through the system and the STM voltage 
is sufficient to find the differential conductance, dl/dV, or the normalized 
differential conductance, d(hiI)/d(hiV) . 



4. Results and discussion 

In this section we show and analyze the current and differential conduc- 
tance as a function of the surface chemical potential fi sur f (ustm = 0). It is 
worth noting that all object atoms are placed directly on the surface and thus 
we assume that their onsite energies are shifted with the surface chemical po- 
tential, Si — > £j + fi sur f- The same is true for the last tip atom and the STM 
electrode. This procedure is justify in our system as the largest potential 
drop is observed for the weakest coupling - in our case this is to parameter 
which is at least 10 times smaller than the other couplings (the tip and the 
substrate are in local equilibrium). All results in this section are obtained 
for the energy unit T = 1. The current and conductance are expressed in the 
units of 2eT/h and 2e 2 /h, respectively, and e.g. for T = 10~ 3 eV the current 
unit corresponds to 10~ 8 A. The other parameters have been chosen in order 



to satisfy the realistic situation in many STM experiments, e.g. [111 . 1 131 [15 

Note that in our calculations we chose fully symmetrical objects in the 
energy space. The first object is a single atom on a surface which is charac- 
terized by the Lorentz-type density of states (DOS). The second one stands 
for a short wire also with symmetrical DOS. It is obvious that in real situa- 
tion such atomic objects on a surface can be characterized by asymmetrical 
DOS e.g. due to different couplings object-surface. In that case the differ- 
ential conductance curves have to be also asymmetrical and it is difficult to 
analyze the role of charged tip in this asymmetry. To avoid this problem we 
concentrate here on objects with symmetrical DOS. 
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4-1. STM tunnelling through a single atom 

For a single atom on the surface the following analytical solutions for V 
and G r matrixes can be obtained: 

• Tip 1, Fig. Eh,: t STM = T 3 ™ , t surf = T sur f and G r = G r n = 

(s-s 1 + ^™ r- f y\ where r STM0 = 2^^1^(8-8,), 

Ysurf _ 2irJ2k(z SU rf l^| 2( K £ ~~ £ k) are energy independent (wide band 
limit approximation). According to the above relations the transmit- 



tance becomes: 

■pSTMO-psnr/ 

\2 , ( r s ™°+r sur f 



T (e) = 1 ; — — (10) 



le - £i 



Tip 2, Fig. [3b (there are two atoms and we assume that the index 
1/2 in all matrixes corresponds to the object/STM tip atom) nonzero 
matrix elements are: ( J 1 stm) 22 = T STM0 , (T S urf)ii = T sur f , and the 
transmittance is expressed as follows: T(e) = Y STM0 T sur ^ \G r l2 \ 2 , where 

\G\ 2 \ 2 = ~, , (11) 



- £i + i^- J [e - Sstm + i^—j ~ *ol 

Tip 3, Fig. [3b, nonzero matrix elements are: (Tstm)ii = V s ™ 1 , 
(fsTM) 2 2 = r 5 ™ and (t surf ) n = T sur f thus the transmittance be- 
comes: 

T(e) = r s ^(r 5 ™ |G^| 2 + r 5TM1 |G^| 2 ) , (12) 

where 

I (e-e 1 + t I^f^l\ f £ I £STM + ^ \ _ ^ 



I (e-e 1 + ( £ _ £STM + z r£^ _ ^ 
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Figure 4: Current-voltage characteristics for the neutral [estm — 0, upper panel) and 
occupied STM atom (estm = — 1 5 lower panel). The broken (thin solid, thick solid) lines 
correspond to tip 1 (tip 2, tip 3) configuration shown in Fig.[3l The other parameters are: 
pSTA/o = rs „r/ = r = 1, V s ™ 1 = 0.1, t Q = 0.1, and e x = /w/, MsTM = 0. 

The knowledge of the transmittance allows us to find the current flowing in 
the system and the differential conductance. 

In order to reveal the role of charged tip we investigate the STM current 
for all described tip models. In Fig. H] the current- voltage characteristics are 
shown for the case of the neutral tip (upper panel) and the occupied one 
(lower panel). The broken lines are obtained for the simplest model where 
the STM tip is represented by a metallic electrode without any tip atom, 
cf. Fig. [3K Thus the broken line is independent on Estm parameter and is 
the same in both panels. The current in this case is monotonic function of 
the voltage. The thin solid (thick solid) lines correspond to tip 2, Fig. [3b 
(tip 3, Fig. [3b). The currents obtained within both tip 1 and tip 3 are 
very similar because the object atom for these cases is coupled directly with 
the STM electrode, cf. the broken and thick lines. In contrast to this, the 
current is not monotonic function for tip 2 and takes nonzero values only for 
Hsurf ~ Estm, cf- thin solid lines, both panels. In the last case electrons can 
tunnel from the STM electrode to the surface only though the tip atom and 
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Figure 5: Differential conductance, dI/dfi SU rf, obtained for the currents from Fig. |4j All 
parameters are the same as in Fig. |4] 

the object atom. This is the reason why the current flows if there are nonzero 
local density of states at both atoms in the window of the chemical potentials. 
This effect is also reflected in the current curve obtained for tip 3 model for 
£stm = — 1 (thick line, lower panel). In this case the current is characterized 
by a local minimum for fi sur f = £stm- The current- volt age characteristics 
obtained for the neutral and occupied tip atom seem to be similar to each 
other, cf. the upper and the lower panel for the broken, thick or thin lines, 
respectively. However, there are very important and crucial for our paper 
differences which become prominent on the differential conductance curves. 

In Fig. [5] we show the differential conductance obtained for the corre- 
sponding current curves from Fig. HI Note that the conductance curves ob- 
tained according to tip 1 are the same at both panels (cf. broken lines) and 
reflect the local density od states of the object atom. As a most prominent 
feature, the conductance turns out to be symmetrical versus fi sur f = for 
the case of neutral STM tip (upper panel) and asymmetrical for charged tip 
atom (lower panel). The explanation of this fact is as follows, i.e. for the 
neutral tip the probability of electron tunnelling from the surface to the STM 
tip or vice versa is the same (the absolute currents for positive and negative 
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Figure 6: Differential conductance, dl j d\x sur f , obtained for tip 3 and for £stm = 0, — 1, —2 
- broken, thin solid and thick solid lines, respectively. All parameters are the same as in 

Kg. a 



voltages are the same, cf. Fig. HJ upper panel). If the tip is charged these 
probabilities are not equal and for electrons it is easier to tunnel from the 
tip to the surface states than from the surface to the STM electrode. 

It is interesting that in the presence of the tip atom we find the negative 
differential conductance, cf. thin solid lines in Fig. obtained for tip 2. 
The negative conductance is related to the discrete structure of the local 
density of states at both (STM and object) atoms (cf. [15( where the negative 
conductance was observed for non-constant surface density of states). If the 
overlap between the tip and the surface states decreases (with increasing 
chemical potential) the negative differential conductance appears, j23|, (see 



24( for other explanation of NDC). This condition can be satisfied only for 



certain range of parameters. For tip 3, thick lines, the negative conductance 
appears only for the case of charged tip (lower panel) and due to the direct 
coupling r 5 ™ 1 , it is not as prominent as for tip 2 (r STM1 = 0, thin lines). 
It is worth to comment on the last result. In real STM experiments it is 
very difficult to observe the negative conductance mainly due to not very 
sharp tip, temperature effects or direct coupling between the STM and the 
surface electrode (not through the object). The last factor gives only constant 
positive value to the conductance (the corresponding current depends linear 
only on the chemical potential difference) and shift the conductance curves 
above the OX axis. Thus in our calculations we omit this coupling and 
concentrate on the charged tip effects. 
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The next intriguing question is whether the energy position of the last tip 
atom, Est M-, influences the asymmetry effect. To investigate this problem we 
chose the third model of the STM system, (viz tip 3, Fig. [3b). In Fig. [6] we 
show the differential conductance curves obtained for Sstm = (the neutral 
tip, broken line), Estm = — 1, and Estm = —2, thin solid and thick solid lines, 
respectively. For Estm < (charged tip) we observe, as before, asymmetrical 
conductance curves versus the zero voltage. Moreover, the more charged 
is the tip the stronger asymmetry is observed. A common feature of all 
curves is the emergence of peaks for /i snr / = Estm and fi sur f = Si, which is 
visible especially for the thick curve, sstm — ~ 2 (the peak for fi sur f = —2 
is related to the energy level of the apex STM atom). This effect could 
be used to test the quality of the STM tip. The apex tip atom leaves its 
fingerprints in the differential conductance curves: First, the conductance 
curves are asymmetrical and second, the conductance peak related to the 
tip atom should appear. Note, that using first-principles calculations, the 
studies of the geometrical, electronic, and dynamic properties of a single 
atoms adsorbed on silicon or tungsten surfaces were reported e.g. in Ref. 
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4-2. STM tunnelling through a single atom - the role of Coulomb interaction 

In this subsection we analyze the influence of Coulomb interaction be- 
tween the STM tip and one-atom object on the electron transport through 
this system, cf. Fig. |3t. The interaction Hamiltonian can be written in the 
following form: 

H int = Uc^CxC% TM CsTM (15) 

Here we restrict our investigation to non-magnetic atoms and do not con- 
sider many-body effects like the Kondo effect, charge-spin separation effect 
and others. It allows us to use a Hartree-Fock like approximation and the in- 
teraction can be captured by renormalizing the on-site energy levels. In this 
case we can substitute E\ by E\ + ^tistm and Estm by Estm + \ n \ where 
n i/STM stands for the charge occupation of the atom on a surface and the 
STM apex atom, respectively. The occupations, uustm, can be obtained 
from the the knowledge of the retarded Green functions according to the 
relation n\ = — f^ rf ImGi 1 (e)de, and similar for Ustm- m this case it is 
impossible to derive analytical solutions for the current or the conductance 
as the retarded Green function depends on both (the surface and the STM) 
atom occupations and is found numerically in the self-consistent way. 
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Figure 7: Differential conductance for the neutral (sstm = 0, upper panel) and occupied 
tip atom {est m = — 2, lower panel). The broken (solid) lines correspond to the Coulomb 
interaction U = to (U = 10io). Tip 3 is considered and the other parameters are the same 
as in Fig. 

In Fig. [7J the influence of Coulomb interaction on the differential conduc- 
tance through one-atom object on a surface is studied. The upper (lower) 
panel corresponds to the neutral, Sstm = (occupied, Sstm = — 2) STM 
tip. For very week interaction, U = t , the differential conductance curves 
are very similar to those ones obtained for U = 0, cf. the broken and thick 
solid lines in Fig. [6] with the broken lines in Fig. [7J For stronger Coulomb 
interaction, U = 10to there are subtle effects in the differential conductance 
curves in comparison with the case for U = 0. The main conductance peak 
related to the surface atom state, for Ustm = 0, is slightly changed for the 
neutral STM tip. For the occupied STM tip the position of this peak is shifted 
towards the negative values of ^stm- It results from the renormalization of 
the on-site atom energy i.e. the STM apex atom is strongly charged which 
shifts the on-site surface atom energy above the fermi level - this leads to 
the conductance peak below fi stm = 0. The results shown in Fig. [7J confirm 
that the Coulomb interactions slightly change the differential conductance 
structure and can be omitted in our studies. 



15 



Moreover, the electron-electron Coulomb on-site interactions are also ne- 
glected in our system as we consider here non-magnetic atoms (e.g. a tung- 
sten _tiPj_an object made of silver, gold or lead atoms and silicon-gold surface, 



cf. ll3|) and, as was mentioned above, do not study many-body effects. 



4-3. STM tunnelling to a short wire 

In order to corroborate the results of the previous subsection we have 
computed the STM differential conductance for a quantum wire on a surface 
(instead of the single atom considered in Sec J4.ip . The wire consists of N = 3 
atoms with the same onsite energies, E\ = 62 = £3, and the nearest neighbor 
couplings, Vi = Vq. The coupling wire-surface is described by the function 
r|J r ^ = T sur ^5ij, and r snr -^ is defined as before. The wire is characterized 
by fully symmetrical local density of states: It is very important in these 
studies - for asymmetrical DOS the conductance is also asymmetrical. In 
our calculations we assume that the STM tip is connected with the first wire 
atom (the results for other connections are similar). The calculated method 
is the same as in the previous subsection thus here we give only the relations 
for the transmittance through three considered systems (cf. Fig. [3] but for 
3-site wire on the surface): 



Tip 1, Fig. Eh: T(e) = r STM0 T sur f (|G^| 2 + \G r 12 \ 2 + \G 



|2> 
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• Tip 2, Fig. [3b (there are four atoms and we assume that indexes: 1,2,3 
describe the wire sites and index 4 corresponds to the tip atom): T(e) = 

T STM0 T surf + | ^ 1 2 + 

• Tip 3, Fig.Eb: T(e) = T STM0 T sur f (\G r u \ 2 + \G r 24 \ 2 + \G r u \ 2 ) + 

T STMl T surf (1^2 + + ^QT^y 

and the retarded Green functions are obtained from the equation of motion 
for these functions. Note, that the transmittance for tip 3 is not a simple sum 
of the transmittance obtained for tip 1 and tip 2 due to Y STMl parameter 
which appears directly in the transmittance relation and also in the above 
Green functions. 

Figure [8] depicts the differential conductance obtained for short quantum 
wire (N = 3 atomic sites) for the neutral (estm — 0, upper panel) and oc- 
cupied tip atom {estm = — 1, lower panel). Such a wire is characterized by 
three-peak density of states which is reflected also in the differential conduc- 
tance curves (cf. the broken line obtained for tip 1, upper panel). The thin 
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-6 -4 -2 2 4 6 
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Figure 8: Differential conductance for the neutral (sstm = 0, upper panel) and occupied 
tip atom (estm = — 1, lower panel). The broken, thin solid and thick solid lines correspond 
to tip 1, tip 2 and tip 3, respectively. Vo = 2, and the other parameters are the same as 
in Fig. [U 



solid (thick solid) lines correspond to tip 2 (tip 3) and these peaks are also 
visible in these cases. For the neutral tip all conductance curves are fully 
symmetrical (upper panel), cf. also Fig. However, for the occupied tip the 
conductance curves are asymmetrical versus [i sur f = (lower panel) which is 
in accordance with the previous results. This asymmetry in the conductance 
is very well visible for the thin solid line obtained for tip 2 (Fig. [HJ, lower 
panel, thin line). Here, for fi S urf = — 3 a local conductance peak appears 
but for Hsurf — 3 there is no such maximum in the differential conductance. 
The conductance curve obtained for tip 3 (thick line) is also asymmetrical. 
It is very interesting and important fact that the occupied tip influences 
all conductance peaks (cf. the thick and thin lines, lower panel) and not 
only the peak near fi sur f = £stm- F° r Msur/ = £stm = —1 one additional 
peak appears which is the same effect as observed for the case of N = 1, 
Fig. This peak is very well visible for tip 2 (thin line) but for tip 3 the 
main conductance peak for fj, sur f = compensates the peak for fi sur f — —I. 
Moreover, as before, the negative differential conductance is observed for the 
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dl/dn surf 




Figure 9: Differential conductance as a function of the surface chemical potential, /i stir /, 
and the position of the last tip atom, £stm for the same parameters as in Fig. [5] obtained 
for tip 3. 



results obtained for tip 2 and tip 3. 

A remaining question is whether the position of e$tm influences in the 
same way all differential conductance peaks. Therefore we plot in Fig. [9] 
the differential conductance as a function of the surface chemical potential, 
fi S urf, and Estm- Here, the symmetry of the conductance for the neutral tip 
is visible for Estm = and asymmetrical shapes of dl/dV for negative Estm- 
It is interesting that the occupied tip influences mainly the conductance for 
negative potentials i.e for fj, sur f < 0. For positive fj, sur f the conductance peak 
is slightly modified. Moreover, in Fig. [9] the evolution of the additional peak 
which appears for fj, sur f = £stm is visible. For small values of Estm this peak 
is compensated by the main conductance peaks (which appear for Estm = 0), 
cf. also Fig. [HJ 

Of course, in the STM experiments it is very difficult (or even impossible) 
to change the tip occupation or Estm parameter. However, it is possible to 
use different kind of tips (blunt or sharp) or other atoms attached to the 
tip. In the last case during the measurement of the same object, one should 
obtain different asymmetries of the conductance as each atom is characterized 
by different Estm parameter. This experimental procedure, however, can 
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destroy the tip or change its geometry which can also lead to asymmetrical 
results. 



5. Comparison with experiment 

In this subsection we compare our theoretical results with the STM exper- 
iment on thin Ag films. Ag was grown on Si(lll)-(6x6)Au substrate held at 
room temperature in an UHV condition. The vacuum chamber was equipped 
with a RHEED (Reflection High Electron Energy Diffraction) apparatus and 
STM (type OMICRON VT). Electrochemically etched and annealed in UHV 
W tip was used. The base pressure was 7 x 10~ n mbar. The average thin 
Ag film thickness was equal to 3 ML of Ag(lll). This procedure yields 
flat Ag islands with (111) orientation, as indicated by the RHEED pattern 
appearance. Figure [TUl upper panel, shows a const ant- current scanning tun- 
neling microscope image of the sample which was used for current-voltage 
spectra measurements. The sample was composed of well separated islands 
with some hexagonal shape. Close inspection of the island's surface mor- 
phology reveals presence of grainy features with diameter of about 1.5 nm. 
Consequently, the surface appears rough, independent of the island height. 
Although RHEED patterns obtained for this sample showed streaky features 
corresponding epitaxially ordered structure, somewhat diffuse width of the 
streaks indicated imperfect order. On Ag islands denoted as A, B, and C in 
Fig. [IUJ upper panel, were measured current- voltage spectra and calculated 



the quantity (dl /dV)/(I /V), which is related to the density of states [271.128 
The thickness of islands determined from the profile lines and counted from 
the wetting layer surface was equal to 1.05 nm for areas A and B and 0.8 nm 
for the island C. Figure [101 lower panel, shows the STS spectra measured on 
Ag islands shown in the upper panel of Fig. [101 Each spectrum in the lower 
panel is the average of 60-90 individual I(V) spectra in order to enhance 
the signal-to-noise ratio. The experiment reported here is a good example 
of a phenomenon frequently occurring during STM measurement practice - 
a sudden variation of the tip length which is a consequence of the tip apex 
rearrangement. Here the scan begins at the bottom, and in the middle the 
average level of the image lowers of 1.5A. Apparently, the tip became shorter. 
In order to keep the tunneling current constant, the STM scanner approached 
the sample. It is also visible that the resolution achieved for the lower part 
is better than for the upper part of the image. Thus the tip switches from 
"sharp" configuration (part I) to the "blunt" one (part II). As the both ar- 
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eas A and B are on the same island, with the same thickness, the electronic 
structure of the sample side of the tunneling junction is identical. We note 
that in following experiment on the same sample we observed also reversible 
variations of the tip length. Figure [10] shows no evidence of touching of the 
sample and one can assume that in course of the scan the tip has lost the 
very top element, presumably a single atom. 

The variation of the topographic image is in accord with electronic struc- 
ture changes shown in Fig. [10J lower panel. The curve A, from the area 
scanned by the sharp tip (part I), is displayed together with the curve B 
measured on the area B with the blunt tip (part II). As a most prominent 
feature we observe two quantum states at the surface voltages V = — 0.4eV 
and V = +0.3eV for both tips (curves A and B) and one additional state for 
V ~ +0.7eV observed only for the sharp tip (curve A). For comparison we 
also measured spectra of the island C with the thickness 1 ML smaller than 
for A and B. It is shown in Fig.[10| lower panel, as curve C. The curves A and 
C are essentially identical. This fact excludes thickness-related changes of the 
thin film electronic structure, namely the Quantum Size Effect (QSE) [29], 
as well as thickness-dependent shift of the Ag surface state binding energy, 
as observed during photoemission experiment 3(|. Although the detailed 



knowledge of the origin of the states observed in STS spectra is not impor- 
tant for the following discussion, we expect that they are relevant to the 
grainy structure of the island. Therefore, one of the possible candidates are 
electronic states of the silver clusters, like the clusters on graphite surface, 



studied with STS in 31 



To analyze theoretically the experimental results we use the method de- 
scribed in the previous section and consider two STM tips: a metallic elec- 
trode, tip 1 (Fig. [3^), and the tip with one apex atom, tip 2 (Fig. [3b). 
The first one corresponds to the STM-tip from part II of our experiment 
whereas the second one is responsible for part I (the tip is sharper). Two 
Ag surface states which are visible in Fig. [TTJJ, lower panel, for V = — 0.4eV 
and V = +0.3eV are generated in our model by means of two coupled atomic 
sites situated on a metallic surface (for our purposes the origin of these states 
is not important). Fig. [TT] depicts the density of states (for two atoms on the 
surface) for the symmetrical (DOS1) and asymmetrical (DOS2) cases which 
are used in our calculations. 

In Fig. [12] the normalized conductance as a function of the surface voltage 
is shown for the symmetrical (upper panel) and asymmetrical (lower panel) 
surface DOS (obtained according to tip 2 model) and for different tip oc- 
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Figure 10: (color online). Upper panel - 30 x 15 nm 2 STM image of an Ag sample deposited 
onto Si(lll)-6 x 6 surface. The rectangles A, B, and C mark areas for which spectroscopic 
/(V) characteristics are analyzed. The scanning parameters are: V s — —1.3V, / = In A. 
Note image contrast difference between part I and part II by sudden STM tip modification. 
Lower panel - normalized dl/dV curves of corresponding areas A, B, and C marked in 
the upper panel. Curves A and B are recorded on the same island, curves A and C are 
recorded on islands with different thickness. Set point before opening the feedback-loop 
was V s = -1.3V, I = InA. 

cupation, Estm = 0, —0.3 and —0.7 - thin solid, solid and thick solid lines, 
respectively. The case of Estm = corresponds to the neutral tip and the 
results are very similar to the differential conductance curves obtained for 
the tip without the apex atom i.e. for tip 1, broken lines. Note, that for 
these cases (metallic tip or the tip with neutral apex atom) the conductance 
curves are characterized by two maxima which reflect only the surface quan- 
tum states. This situation corresponds to the blunt STM tip used in our 
experiment (part II). The thick solid lines in Fig. [12] correspond to the case 
of occupied tip, Estm = —0.7, and the additional peak for V = +0.7eV is 
visible. (Remark: the STM tip states which are below the Fermi level, on 
the conductance graphs are visible for positive surface voltage). It is high 
probably that the origin of this peak is related to the apex atom energy level 
which is below the Fermi energy and is occupied (the surface states in the 
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Figure 11: Density of states obtained for two coupled atoms, N = 2, on a surface. The 
broken line corresponds to z\ = £2 = —0.05, Vq — 0.35 (symmetric case), the solid line is 
obtained for £1 = —0.35, £2 = 0.2, Vq = 0.2 (asymmetric case). 



experiment are the same for both parts but the tip has changed). The similar 
effect appears in Fig. [5] or Fig. [6j where for charged tip additional structure 
in the conductance was observed. Thus, the broken curves in Fig. [T2] rep- 
resent the STM tip from part II and the thick solid lines correspond to the 
tip from part I. Moreover, for asymmetrical DOS, lower panel, the qualita- 
tive comparison of the differential conductance with the experiment is better 
which suggests that the intensity of Ag surface states (below and above the 
Fermi level) are different, cf. Fig. [TT], DOS2. 

It is interesting that the normalized differential conductance peaks are 
not proportional to the surface DOS and the results obtained for symmetri- 
cal and asymmetrical DOS do not reflect strong asymmetry which appears 
for DOS2, cf. Fig. [TTJ It means that the intensity of d(\n I) / d(ln V) peaks are 
not directly proportional to the surface DOS (the positions of these peaks are 
related to the surface or STM DOS), cf. also (27| . Only for nearly constant 
tip density of states the normalized differential conductance is proportional 



to the sample DOS, [28(. Thus, using STS results it is difficult to distinguish 
between the surface and the STM tip states due to their energy convolution. 
One possible way is to scan two or more different areas of the surface (char- 
acterized by different structure or different atoms) - if there are some states 
which appear for both areas (for the same voltage) there are probably the 
STM states. 

To conclude, the STM states influence the conductance behavior and dis- 
turb the intensity of the conductance peaks. For charged STM tip, Estm < 0, 
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Figure 12: Normalized differential conductance d(ln/)/d(ln V) versus the surface voltage, 
V = Ustm — Hsurf, obtained for symmetrical DOS (upper panel) and asymmetrical DOS 
(lower panel) and for different apex atom energy (according to tip 2 model): sstm = 
0, —0.3 and —0.7 - thin solid, solid and thick solid lines, respectively. Surface DOS are 
shown in Fig. QT] The broken lines correspond to the metallic tip without the apex atom 
(tip 1 model). The other parameters are: r 5TM0 = 0.5, r STM1 = 0, t = 0.1, and 

USTM = 0. 



it leads to the symmetry braking in the differential conductance, dl/dV, 
which is also visible on the normalized differential conductance curves, cf (In /) / d(ln V) . 
It is worth noting that the position of the additional peak on the differential 
conductance curve, Fig. [T2J obtained for V = +0.7eV is in good agreement 
with the results reported in Ref. Q, Fig. 6, where the authors investigated 
Au(lll) surface by means of the tungsten tip and observed a special peak for 
V = +0.6eV which, in their opinion, was a halo of the tip structure. In our 
studies the position of the tip-induced peak is very similar which indicate 
that in experiments with sharp tungsten tip a local maximum around the 
sample voltage 0.6 — 0.7eV should be observed. This is also consistent with 



the STS results on Cu(lll) surface, [32|, where the differential conductance 



peak for the sample voltage V = 0.5 ± 0.3eV appears and is identified as a 
halo of the tungsten sharp STM tip (fee pyramid configuration, cf. also 261]). 
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6. Conclusions 



In summary, using the Green's function technique for the tight-binding 
Hamiltonian the STM transport properties through surface states have been 
studied. It was shown that sharp STM tip, represented by a kind of bending 
wire, is extra occupied. In this case the analytical relation for the Green's 
function needed to obtain the local charge has been obtained for arbitrary 
wire length, N, Eq. [2j 

The current and the differential conductance have been studied for three 
models of the STM tip and for two objects on the surface (single atom and 
short wire). The analytical formulas for the transmittance and the retarded 
Green functions have been obtained. For the neutral STM tip the conduc- 
tance has turned out to be symmetrical versus the surface chemical potential, 
(cf. Fig.[5j upper panel) and asymmetrical for charged tip atom (Fig. [5J, lower 
panel and Fig. [9]). Thus for occupied tip all differential conductance curves 
are asymmetrical, although the density of states of investigated objects are 
fully symmetrical in the energy space. Additional peak in the conductance 
has been observed for the STM voltage which corresponds to the onsite en- 
ergy of the last tip atom. It means that this peak is a feature of the STM 
tip and thus it should be observed for different objects at the surface. 

To confirm our theoretical calculations the STM experiment on Ag islands 
with two different tips has been carried out. Additional peak in the conduc- 
tance curve for +0.7eV has been observed for sharper tip which is a halo of 
the tip quantum state. This peak has not been registered for the blunt tip, cf. 
Fig. HOI Our theoretical calculations obtained for charged STM tip, Fig. [121 
are in good agreement with the experimental results. Alternatively, such an 
experiment may be also performed with metallic molecules or clusters on a 



surface, e.g. [33, |34|. It has been also shown that the normalized differ- 
ential conductance peaks are not directly proportional to the surface DOS. 
Moreover, the results obtained for symmetrical and asymmetrical DOS do 
not reflect asymmetry in the peak intensity of the surface DOS. 
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